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ABSTRACT 

We investigate the potential of the Square Kilometer Array Telescope (SKA) to con- 
strain the sound speed of dark energy. The Integrated Sachs Wolfe (ISW) effect re- 
sults in a significant power spectrum signal when CMB temperature anisotropics are 
cross-correlated with galaxies detectable with the SKA in HI. We consider using this 
measurement, the autocorrelation of HI galaxies and the CMB temperature power 
spectrum to derive constraints on the sound speed. We study the contributions to the 
cross-correlation signal made by galaxies at different redshifts and use redshift tomog- 
raphy to improve the signal-to-noise. We use a analysis to estimate the significance 
of detecting a sound speed different from that expected in quintessence models, find- 
ing that there is potential to distinguish very low sound speeds from the quintessence 
value. 

Key words: cosmological parmeters - large-scale structure of the universe - cosmic 
microwave background - radio lines: galaxies 
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1 INTRODUCTION 

Current observations of the Cosmic Microwave B ackground 
(CMB) anisotropy spectra jSpergel et al.l |2006| . hereafter 
WMA P3) combined with large scal e structure ( LSS) sur- 
veys (|Percival et al.l I2OOII : Tcgmark 'et all |2004|) and su- 
pernova Type la (SNIa) experiments ( Riess et al.l 1 19981 : 
iPerlmutter et af][l999l ) have provided strong evidence that 
we are living in an Universe dominated by an unknown and 
exotic form of matter which has been called dark energy 
(DE). 

Further evidence for this dark component is pro- 
vided by the detection of a positive cross-correlation 
between the CMB temperature anisotropy and several 
LSS su rveys through the Integ r ated S achs Wolfe (ISW) 
effect jAfshordi. Loh fc Strauss! [20oi: iNolta et all 12004 



scalar field models cannot be a priori excluded. In fact, 
the small (but non-vanishing) value observed for the 
dark energy density (pde) allo ws the possibility of study- 
ing d ynamical models (e.g. Zlatev. Wang fc Steinhardtl 
19991: lArmendariz-Picon. Mukhanov fc Steinhardtl I2OO0I : 



Amendolall2000l ) 



Boughn fc CrittendenI l2004l : iFosalba fc Gaztanagal I2OO4 



Cabre et al. 20061 ). This effect arises when CMB photons 



streaming across the dark energy dominated Universe inter- 
act with the decaying gravitational potential wells associ- 
ated with the foreground LSS. This interaction imprints a 
secondary anisotropy on the primordial CMB spectrum at 
large scales and it can be used to probe the nature of dark 
energy. 

Current models of dark energy include a cosmolog- 
ical constant (A), with constant energy density filling 
space homogeneously, and dynamical scalar field mod- 
els, such as quintessence, whose energy density varies 
with time and space. Although the current concordance 
model, involving A and cold dark matter in an FRW 
Universe (ACDM), is consistent with all observations. 



A common way of describing DE models is through its 
equation of state parameter, w = pde/ Pde, where pde is 
the pressure and pD e is the energy density. The cosmological 
constant model corresponds specifically to u; = — 1, whereas 
a value of w 7^ —1 is associated with dynamical models. 
The equation of state parameter affects both the evolution 
of density fluctuations and the background expansion of the 
Universe. 

Another characteristic that defines the nature of a gen- 
eral matter component is its effective sound speed, Ci. For 
a perfect fiuid with adiabatic perturbations there is a sim- 
ple relation between the background and the sound speed: 
Ca = p/p. In imperfect fiuids, like scalar field models of DE, 
however, dissipative processes generate entropic perturba- 
tions, and the sound spee d is given by the more general 



IS given 

relation: ci — Spde/Spde (|Hulll99j ). The sound speed of 
DE energy thus affects the growth of density perturbations. 

Both the equation of state parameter and the sound 
speed determine the clustering properties of the DE fluid. 
Single scalar field models, or quintessence, are characterised 
by having a speed of sound very close to the speed of light, 
and DE clustering is determined by w alone, occu rring only 
on sc ales equal to or larger than the horizon size (|Ma et al.l 
I1999I ). For alternative sound speed models, with < 1 (in 
units of the speed of light, c), DE clusters on sub-horizon 
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scales which are determined by the value of Cs- An exam- 
ple is the scalar field model named k-essence, characterised 
by having a non-stan dard form of kinetic ene rgy with an 
evolving sound speed (|Armendariz-Picon et aLi .2000). 

The ISW effect presents itself as a probe of the 

properties of DE (|Corasaniti et all 120031 ) and several 
groups have provided new ways of constraining the 
DE parameters using the cross-correlati on between the 
CMB and different LSS surveys (e.g lAfshordU |2004 
Hu fc Scraiitoiil l2004l: iBean fc Dorl l2004l: IPoeosian et all 



2OO5I : Corasaniti. Giannantonio fc Melchiorrill2005l ). In this 
paper we present the prospects of constraining the proper- 
ties of DE by using an Hl-galaxy survey that may be attain- 
able with the future Square Kilometre Array (SKA) radio 
telescope using the ISW effect in cross-correlation and the 
Hl-galajcy autocorrelation. In particular, we investigate the 
prospects for detection of a DE speed of sound different to 
that of quitessence (i. e., ~ !)■ 

The outline of the paper is as follows; in ij2] we give 
an overview of the theory of the ISW effect describing the 
effects of the DE parameters on the ISW-galaxy and galaxy- 
galaxy power spectra. In i|3] we introduce and discuss the 
properties of a model SKA HI redshift distribution and the 
impact on ISW experiments. In section ^we deal with the 
numerical implementation of the theory, the signal-to-noise 
calculations and the forecast for model detection. Finally, 
the results are presented and discussed in ^ and they are 
summarised in il6l 



where A^(fc) = Ps{k)/2'K'^ is the matter power spec- 
trum per logarithmic interval of wave number (k) today and 
Ps{k) = {\5{k)f). The filter functions /f (fc) and /f (fc) in 
Eq. (Is} correspond to the Hl-survey and ISW spectrum re- 
spectively. The amount of cross-correlation signal will de- 
pend on both the distribution of the galaxy selection func- 
tion and the change of the gravitational potential with time. 
This dependency is manifested through the two filter func- 
tions which are defined as integrals over redshift weighted 
by spherical Bessel functions, ji{kr{z)), as 



fl{k) 



dN 

bHi{z)——D(z) jt{kr{z)) dz , 
dz 



c^k^ 



dz 



ji{kr{z))dz, 



(4) 
(5) 



where Ho is the value of the Hubble constant today and Qm 
is the matter density today in units of critical. The function 
g{z) is the growth suppression factor which is related to the 
linear growth of matter perturbations as g{z) = (H- z)D(z), 
with D{z) given by 5m{k, z) — D{z)Smik, z = 0). 

Similar to Eq. the galaxy autocorrelation coeffi- 
cients (the galaxy angular power spectrum) can be obtained 



= 4tv 



jt (kr) 



A^(fc) 



(6) 



2 THEORY 

2.1 CMB— LSS cross-correlation and galaxy 
autocorrelation 

Fluctuations in the matter distribution of the large-scale 
structure can be expressed in terms of the projected frac- 
tional source count of the mass tracer 



5N 



bHi(z) '^5m{z,h) dz, 



(1) 



where bHi is the linear bias parameter of an Hl-galaxy pop- 
ulation, 5m is the matter density contrast (= Sp/p), and 
dN /dz is the normalised redshift distribution of HI galax- 
ies. 

On the other hand, the temperature fluctuations arising 
from the ISW effect are expressed as the change in gravita- 
tional potential over conformal time (or comoving distance) 



ST r-^" 

— (n) = -2 / dr^{rn) , 

JO Jo 



(2) 



where $ is the Newtonian gravitational potential, r is the 
comoving distance running from today to the epoch of de- 
coupling and the dot represents derivatives with respect to 
r. 

The amplitude of the galaxy-temperature cross- 
correlation for each mode in Fourier space is thus given as 



An 



dk 



f?{k)fJ{k)Al(k), 



(3) 



2.2 The effect of DE cosmological parameters 

The background expansion of the Universe, determined by 
the Hubble parameter H{z), is affected by the density of 
the matter components and by w. The comoving distance 
appearing in the previous expressions is defined as 



r = / dz'/H{z'), 
Jo 



(7) 



where, for a fiat Universe and in an epoch where radiation 
is negligible, H{z) is given by 



Ho 



= \ flmil + z)'^ + Qde exp 



3 / dz 





l + w{z) 



1 + z 



(8) 



where we allow for a time-dependent w. For a constant equa- 
tion of state parameter, Eq. ((Sjl simplifies 



H{z) = Ho \ (1 + z)3 + QoE (1 + 2)3(1+"') 



(9) 



The dark energy term in Eq. (|9]) has also an effect on 
the growth of matter fluctuations. For quintessence mod- 
els, on scales smaller than the horizon, dark energy can be 
con sidered smoot h and the matter overdensity grows as (see 
e.g. |Peebles|[l980l ) 

3 



+ 2HS 

rr 



-H\l - Q.K)5rn = 0. 



(10) 



Eq. (|10p is usually rewritten as a function of the linear 
growth factor, D{z,w), and can be solved numerically to 
compute, for example, the power spectra in the above equa- 
tions for each w model. 

More generally, the introduction of DE perturbations, 
characterised by Cg, also has an impact on the evolution 
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of CDM overdensities. This effect is, however, less signifi- 
cant for smaller w and it vanishes as w —1 (in fact, it 
is only in this limit when one could safel y neglect DE per- 
turba tions). In the synchronous gauge (Ma fc Bertschingeij 
1 19951 ). the evolution of the density and velocity perturbation 
of a general matter component (denoted by subscript i) in 
the frame como ving with the DE (den oted by the circumflex 
(")) is given by (jWeller fc Lewisll2003l ) 

Si + SHi&l, - + 3^(1 + w^)v^/k) + 

(1 + Wi)kvi = —3(1 + Wi)h' 



(11) 
(12) 

where H is the conformal Hubble parameter, Vi is the 
veloc ity and h' is the syn chronous met ric perturba- 
tion ( Ma fc Bertschingerl Il995l ) . As shown in iBean fc Dord 
the effect of increasing the sound speed of the dark 
energy is to reduce the growth of structure resulting in an 
increased ISW signal on the largest scales. On the other 
hand, as Cg decreases, CDM perturbations larger than the 
DE sound horizon collapse, and the ISW effect is reduced. 

Eqs. (Ilip and (|12p show how density perturbations are 
fully characterised by both w and c^. The solution to these 
equations can be obtained numerically, provided that we 
know these two parameters. In this work, we w i ll con sider 
these to be constants, as justified in iHannestadl l|2005h and 
references therein, although the formalism can be easily 
adapted for evolving models of dark energy. 



3 THE SKA HI-GALAXY REDSHIFT SURVEY 

3.1 The SKA for ISW experiments 

Direct observation of HI emission from neutral atomic 
hydrogerQ at high redshift is limited by the poor sen- 
sitivity of current radio telescopes. Direct detections 
have been only achieved up to red shift z ~ 0.18 
(IZiwaan. van Dokkum fc VerheiienI l200lh . and HI gas at 
higher redshift has only been probed through the absorp- 
tion in the quasar optical spectra by damped-Lya (DLA) 
systems, in which the bulk of the neutral gas i n the Universe 
at hi gh redshift is believed to be contained jPeroux et al.l 
I2OO3I ). 

The SKA will have an unprecedented sensitivity pro- 
vided by its large collecting well as a large (~ 10 
deg^) instantaneous field of view (FOV). The telescope's 
high sensitivity will allow the detection of HI galaxies lo- 
cated at redshifts that include the beginning of DE domina- 
tion era, z ~ 0.5 - 1, and beyond. This will provide an oppor- 
tunity to map the sky in different redshifts slices (redshift 
tomography), and therefore to constrain more accurately the 
value and evolution of the DE parameters. In addition to a 
fast surveying capability, an extended FOV will also pro- 
vide a very wide survey area (almost the entire sky) which 
is needed to probe the largest scales required by an ISW 
experiment. All the above will result in the detection of a 
large number of sources which will reduce statistic errors 
greatly, and will provide measurements of cross-correlation 
and galaxy autocorrelation with large signal-to-noise. 



^ 21 cm line, produced by a change from a parallel to antiparallel 
spin configuration for the proton and electron. 



3.2 The Hl-galaxy survey 

Our current knowledge of the expected HI galaxy redshift 
distribution is very li mited. A detailed analysis in this re- 
gard is presented in lAbdalla fc Rawlingj l|2005l . hereafter 
AR2005) where the HI mass function, dn/dlogioMm, is in- 
tegrated over HI mass 



dN_ 

dz 



dn 



dV 



dlogioMjii dz 



dlogioMm , 



(13) 



where V is the comoving volume and Mm{z) is the limiting 
mass that can be detected by the SKA under certain as- 
sumptions. The HI mass function inside this integral is very 
likely to evolve with redshift due to the fact that the content 
of neutral Hydrogen is affected by many astrophysical pro- 
cess such as star formation, supernova explosions, etc. The 
models for the evolution of this HI mass function provided 
in AR2005 assume the Schechter paramctrisation measured 
by the HIPASS team at redshift z ~ (Zw aan ct al. 2003) 



dn 



exp 



(14) 



d(AfHi/M*i) 

with a = —1.3. The value of the normalisation, 8*, is set 
according to the total density of neutral gas Qhi found 
at redshift z, whic h is obtained from DLA observations 
(jPeroux et al.ll2003l ). The break of the mass function, Mhj, 
is, however, dependant on the model for galaxy formation 
and the evolution of the HI mass. In AR2005, three models 
are presented for this break and we adopt their preferred 
and more realistic Model 'C, which takes star formation 
into account. 

The limiting HI mass in Eq. (|13|l achievable with the 
SKA is dependant on several observational parameters. For 
a survey where all redshifts are accessible through single 
pointing, a FOV frequency dependence of f'^, a detection 
threshold of S/N = 10, and in 4-hour integration time, we 
use the redshift distribution expressed by the following fit- 
ting formula 



dN 
dz 



■ A z exp 



jz - z,Y 
2Az^ 



(15) 



with A = 2.52 X 10^ Zc = 0.211 and Az = 0.461. The units 
are galaxies per square degree. Fig. [1] shows the resulting HI 
galaxy distribution. 

We define our normalised galaxy selection function as 
0(2;) = bm{z) X dN /dz, where the tilde represents the nor- 
malised distribution over the redshift of interest. The bias 
factor used for our analysis (depicted in the lower panel 
of Fig. [T]) i s adap ted fr om the pa rametrised halo model of 
IHu fc J aini (|2004l ) (as in lHu fc Scra nton 2004), and modified 
at low redshift to match recent HIPASS clusteri ng measure- 
ments (|Basilakos et al.ll2006l : iMever et al.ll2007l ). We varied 
the bias model and found it has little effect on our estimate 
for the significance of sound speed detection. 



4 CALCULATIONAL METHOD 
4.1 Numerical implementation 

To test the sensitivity of the cross-correlation to vari- 
ations in the DE parameters we compute numerically 
the relevant power spectra (Eqs. ((3} & (|6]) plus CMB 
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Figure 1. Top panel: The number of Hl-emitting galaxies per 
square degree as a function of redshift expected to be detectable 
with the SKA for an integ ration time of 4 hours. See text and 
lAbdalla fc Rawlingj 1I2OO5I ) for details. Also in grey are repre- 
sented the redshift slices for S/N measurements (see lUJ. Bot- 
tom panel: Our mod el for the bias parameter adapted from 
IHu fc ScrantonI |20oi) (dashe d) to match recent HIP ASS clustcr- 
ing m easurements (solid line) l lBasilakos et al. I I2OO6I : IMever et al.l 
l2007t ). 



power sectrum) using a modified version of CAM^l code 
(|Lewis. Challinor fc Lasenbvl l2000l ). The code takes our 
galaxy selection function combined with the matter over- 
density, 5m{k,z), and solves the line-of-sight integral of 
the resulting ISW source term over conformal time (as in 
Eq. (jUlfl). The output contains the coefficients for the CMB 
auto-power, Cf^, the cross-power with galaxies, Cf"^ and 
the galaxy auto-power, Cf^. 

When considering different models of DE we want to 
isolate the large-scale variations of the CMB spectrum (i.e., 
the contribution from the ISW alone). For this reason we 
fix the angular diameter distance to recombinatiorQ and 
the physical matter density, Qrnh^, to be consistent with 
WMAP3 best fit. To create a family of models, we choose 
values for w and h{= Ho/ WO) that keep the angular diame- 
ter distance fixed, and then vary the DE sound speed, Cg. To 
generate the CMB data, we set the optical depth to reioni- 
sation, r — 0.09, the primordial power spectrum amplitude, 

= 2.02 X 10^', the scalar spectral index, Us — 0.951, 
and the ratio of the initial tensor and scalar power spec- 
trum amplitudes, T/S = 0.55 {WMAP3). 

In order to increase our signal-to-noise, we divide the 
selection function into several redshift slices (bins) and con- 
sider the cross-correlation in each bin as an independent 
measurement. This assumes that the SKA will be able to 
measure the redshift of each unresolved galaxy accurately 



^ |http : //camb ■ Inf o| March 2006 version. 

^ Or similar, as CAMB computes 5m, {k, z) from the evolution equa- 
tions (Eqs. jlip and l|12p ) and considers instead the primordial 
curvature power spectrum in Eq. 

* This is well constrained by the position of the first acoustic 
peak in the WMAP observations. 



through the HI line. Within this approximation, we can 
then add up the values of (S/N)^ for each redshift bin in 
order to obtain the total signal-to-noise of the experiment. 
The width of the bins is chosen carefully. Wide bins pro- 
duce larger cross-correlation signal through the ISW filter 
(Eq. (O), as there is a greater variation of the potential. 
This, however, increases the total variance as there are fewer 
measurements along the redshift space. Narrow bins, on the 
other hand, suffer from a slow-varying gravitational poten- 
tial and shot noise. Considering the above effects, we choose 
in our calculations an optimal bin size of Az — 0.2 (see 

Fig. DO). 

4.2 Signal-to-noise measurement 

The ISW-galaxy cross-correlation signal is affected by var- 
ious sources of noise such as the cosmic variance, Poisson 
noise of t he galaxy surve y, CMB pixel noise and redshift 
errors (see lAfshordll2004l for more details). 

The signal-to-noise ratio squared for a cross-correlation 
measurement at multipole I due to a redshift bin, i, is addi- 
tive 



{S/Nf = J2 



(16) 



where the sum extends over redshift bins and multipoles. 
The expec ted error ACf can b e obtained from basic prin- 
ciples (e.g. lAfshordi et al]l2004 ) 

1 



AC; 



9T 



(2^+i)./.i; 



^i(c7f )^+crcr 



(17) 



where f^^^ is the smallest of the CMB and galaxy survey sky 
coverage fractions and the coefficients Cf^ and Cf^ include 
the noise power spectra. For the galaxy field, the source of 
noise comes from the Poisson fluctuations in the number 
density 



CI' = Cf" + \ln\ 



(18) 



where h\ is the galaxy number per steradian in the redshift 
bin of interest. For the CMB, the noise contribution to the 
temperature depends on the beam window function and the 
pixel noise of the experiment (see below) . 

4.3 Significance of detection 

In order to determine the signiflcance level of a DE sound 
speed detection away from the reference model = 1 for a 
given ui, we use the statistic 

-2iogr^/ = E'''^'^T^^^"''^'^"^'lr'^^ (19) 



[AC/(u>;ci = l)]' 



where L is the likelihood of observing the power spectra for 
a given set of cosmological parameters. In this analysis we 
assume that, by the time of SKA construction and opera- 
tion, future independent experiments will provide accurate 
measurements of parameters other than c^. We will explore 
potential degeneracies and the possible evolution of in a 
follow-up work. 

The superscript A in Eq. (|19p refers to the different 



observables that constitute our data set, i.e., Ci , C?® and 
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Table 1. CMB Experimental Specifications. Detector sensitivity 
aT is in per FWHM beam and the weight given to one channel 
is wx = {^^bcamO'/TcMB ■ the CMB map we take a sky 



0.1 



coverage of /sky = 0.8 and thus set 



ky 



0.8. 



Frequency (GHz) 


fbeam (arcmin) cry (^K) 


s 


143 


8.0 6.0 


u 


217 


5.5 11.4 











Cf The errors for the CMB temperature anisotropies and 
a full-sky {f^^^ = 1) galaxy autocorrelation measurement are 
respectively 

2 



AC 



(2f +l)/sky 



(2£+l) 



icf) 



(20) 
(21) 



99T 

^max 



We extend the sum in Eq. (I19|) up to t , 
for the cross-correlation and autocorrelation experiment^, 
and set £^ax = 2000 for the CMB, although the exact value 
in this case is not critical for sound speed separation through 
the ISW. Note that we consider independent redshifts mea- 
surements and therefore do not include the cross-correlation 
between bins in Eq. (|19p . 

The CMB temperature power spectra, Cf^, include the 
contribution to the pixel noise of the CMB experiment. The 
observed spectra, i.e. signal plus noise, are 



Ci ~ Ci + w. 



1 D-2 
T ■ 



(22) 



where Bj = exp(-^(^ + l)eLam 



£ = expi^— {.(^t -r j-jpjj^j^j^/ 8 In 2) is the window func- 
tion of the (assumed) Gaussian beam, and Sbeam is the beam 
FWHM. The factor wt accounts for the detector noise in the 
temperature measurements. We assume a CMB experiment 
that measures only temperature anisotropies in two high fre- 
quency bands, 143 and 217 GHz (Rocha et al. 2004 ). In this 
case, the value of wtBJ in Eq. (|22|l is the sum of wtBJ 
for each frequency channel. The paramete r s for t his type of 
experiment are adapted from lRocha et ahl l|2004h and listed 
in Table [U 

Using a similar approach to lPogosian et ah! |200^, we 
combine our three experiments in order to measure the total 
significance of separation of sound speed models, (S/N)'^, 
using Eq. ^ 



(S/N) 



+ X 



(23) 



5 RESULTS 

5.1 The SKA-CMB cross-power spectrum: 

Firstly, we show the effect of the sound speed and the equa- 
tion of state parameter on the cross-correlation spectrum of 
the SKA Hl-galaxy field with the CMB. This is illustrated 
in Fig. [2] where we consider different values of w and and 
a single selection function which covers the redshift range 
from z = to z = 2. Our experiments have shown than 



^ In this analysis we ignore CMB polarisation data. 

^ Beyond I ~ 300, the galaxy autocorrelation becomes dominated 

by the Poisson term. 
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Figure 2. Power spectrum from cross-correlation of galaxy over- 
densities and CMB anisotropies for a single redshift bin < ^ < 2. 
For the different models, the angular diameter distance to CMB 
is fixed as well as Qrn.h'^ • The cross-correlation becomes more sen- 
sitive to the sound speed as w increases. 



this redshift configuration does not provide the optimal S/N 
which is needed for model comparison. This diagram, how- 
ever, helps to illustrate the effect of w on the prospects of 
distinguishing a 7^ 1 scenario. 

Fig. [3] shows the cross-power spectrum of each red- 
shift bin referred to in Fig. [T] with the CMB for a ref- 
erence w = —0.8 and quintessence model. As expected, 
the cross- correlation signal at higher redshifts peaks at 
smaller angular scales (the galaxy correlation length looks 
smaller for deeper objects). Because these scales are least 
affected by cosmic variance, we expect high-redshift objects 
to contribute greatly in model separation through the cross- 
correlation. This is seen in the lower panel of the figure where 
we show the change in the ratio of the spectra for = 0.01 
and Cs = 1 models as a function of redshift. The angular 
scales that provide the largest model separation (the min- 
ima in the curves) get smaller (shift to higher ^'s) as we 
move to higher redshifts. 

5.2 The galaxy autocorrelation: C*® 

We obtain additional information when we look at the an- 
gular autocorrelation of the HI galaxies. This is presented 
in the top panel of Fig. |4] as a function of redshift for 
the quintessence model. This model is compared with the 
Cg — lO"** scenario in the lower panel. On small angular 
scales, the value of the sound speed has little impact on the 
galaxy perturbations and we recover the smooth scenario of 
DE. This effect reduces our ability for model separation us- 
ing the autocorrelation on the smallest scales, although it 
is slightly compensated by the shift of the curves towards 
larger I's caused by projection effects (as mentioned ear- 
lier). The suppression of the curves' amplitude with redshift, 
caused by the onset of DE domination occurring only at re- 
cent times, affects also our ability for model separation. We 
include, in addition, the ratio of the = to the smooth 
case for the lowest redshift bin, and find the discriminating 
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Figure 3. Top panel: The cross-power of the binned HI galaxy 
field with the CMB for a reference model ui = —0.8 and = 1. 
Bottom panel: Ratio of the = 0.01 and = 1 spectra. In 
both panels, the thick solid lines represent the closest and farthest 
redshift bins, and the intermediate bins are represented by dashed 
Unes. 
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Figure 4. Top panel: The angular autocorrelation of the binned 
HI galaxy field for a reference model w = —0.8 and = 1. 
Bottom panel: Ratio of the = 10"** and = 1 spectra. Bins 
are represented as in Fig. |3] The additional grey line at the top 
of the panel gives the ratio of the = and = 1 spectra. 
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Figure 5. (a) For the ISW effect in cross-correlation: The dilfer- 
ential (S/N)'^ per unit redshift for various DE sound speeds as a 
function of redshift. (b) The total (S/N)'^, resulting after adding 
up the contribution of each redshift bin with z sC Zmax • 



we show the (S/N) for a /f^^ = 0.8 survey in a «; = —0.8 
cosmology and different values of the sound speed as a func- 
tion of redshift. Fig. [5^ gives the differential (S/N)^ per unit 
redshift, whereas Fig. [Sja shows the total (summed) (S/N)^ 
up to redshift Zmax- Although the cross-correlation signal-to- 
noise peaks at the epochs of DE domination, the results on 
Fig. [SJd suggests that higher redshifts are required in order 
to attain a reasonable (S/N) for ^ models. 



power to be almost scale-independent as, in this limit, DE 
clusters on all scales. This is an advantage that will become 
visible when we present the significance for detecting Cs ^ 
models. 



5.3 Cross-correlation signal-to-noise 

Next, we want to investigate the effect of the sound speed 
and redshift on the cross-correlation signal-to- noise. In Fig. [5] 



5.4 Significance of DE sound speed detection 

Finally, we present our main results regarding the capabil- 
ity of the SKA to separate different sound speed models of 
DE. The significance of separation between different models 
of sound speed with respect to the quintessence model is 
shown in Fig.[6]for aw = —0.8 model and all the other cos- 
mological parameters kept fixed. The significance increases 
gradually as the sound speed approaches = 0, but falls 
off rapidly as 1. In the figure we show the effect of 

combining the three separate experiments we referred to 
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logio(c^) 

Figure 6. The significance of separation between cosmologies 
with ^ 1 and quintessence with all other parameters fixed and 
fiducial value w = —0.8. Information of the galaxy autocorrela- 
tion combined with the ISW— galaxy cross-correlation increases 
the S/N as — > 0. Sound speed measurements are improved by 
the use of higher redshift (z > 1) objects. 

in § 14.31 The CMB experiment gives poor constraints on 
its own, whereas a detection of Cs = 0.01 at S/N = 2 is 
achieved by cross-correlation using the full redshift survey 
and a,t S/N = 3 when the galaxy information is included. 
As Cs 0, the effect of adding the galaxy autocorrelation 
becomes important, reaching S/N = 5.5 with all the exper- 
iments combined. Although the ISW effect arises mostly at 
low redshifts, deeper measurements can improve the signal- 
to-noise in cross-correlation by exploiting the contributions 
from higher-^ multipoles (see Fig. [3| . 

The significance of detection of a zero sound speed 
model, Cs — 0, against a quintessence model is depicted in 
Fig. [7] for different given values of the equation of state pa- 
rameter, w. Within the current uncertainty in the value of w 
(within 10 per cent of a ACDM model) , we find that a detec- 
tion of Cs = will have a low significance of S/N 1 for the 
w = —0.95 case, but this situation improves at S/N — 2.4, 
for AW — —0.9 Universe. These results give a first estimate of 
the significance 5'/A^ but a full statistical analysis is required 
before strong conclusions can be drawn. We also note that 
although current observations do not favour models with 
w > —0.9, constrai nts could be relaxed by cons idering alter- 
native results (e.g. iDe La Macorra et al]|2006l ). 

For completeness. Fig. [5] summarises our previous re- 
sults by exploring the significance S/N boundaries for sound 
speed separation for given values of w. As expected (see 
Fig. [21), the significance increases as we move along models 
with decreasing sound speed for increasing w. 



6 CONCLUSIONS 

In this paper we have studied the potential of the SKA to 
constrain different models of dark energy and, in particular, 
its speed of sound by combining CMB data and an Hl-galaxy 
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Figure 7. The significance of detection of a DE sound speed 
(?s = 0, for a quintessence model with all other parameters fixed, 
for a given w. 




Figure 8. Significance boundaries for detection of sound speed 
for given w, in a quintessence model Universe. Models along the 
boundary lines are detected with the marked significance. The 
S/N increases as we move along models of increasing w and/or 
decreasing c^. Models with w —0.9 can be only detected as 
— > (see text). 

survey attainable with the SKA. For a constant equation of 
state parameter, we find that the power spectrum of the 
cross-correlation is only sensitive to the sound speed as w 
separates from the cosmological constant limit (w = — 1). 
For a reference value of w = —0.8, low redshift surveys 
(z ~ 0.5) only have the potential for sound speed model 
separation at the largest scales where cosmic variance dom- 
inates, and we find that higher redshift objects (z 2) are 
still needed to increase the total signal-to-noise in the cross- 
correlation. This is particularly important as Cs — + where 
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the total cross-correlation S/N is the lowest {S/N 5). Our 
results have not shown a strong dependence on the chosen 
HI bias model. 

In order to forecast the significance of sound speed sep- 
aration, we have provided a simple analysis containing 
information from three separate measurements: a CMB tem- 
perature experiment, the Hl-galaxy autocorrelation mea- 
surement and the ISW-galaxy cross-correlation. We find 
that we can only constrain variations of sound speed that are 
orders of magnitude away from the reference value, ~ 1- 
For very small sound speeds, useful constraints can be ob- 
tained from the galaxy autocorrelation measurement. As 
Cs — > 0, models with w within 10 per cent of ACDM can 
be separated with at a significance level 2 < S/N < 3. 

A Fisher matrix analysis including temperature, po- 
larisation and galaxy spectra combined together in a sin- 
gle data covariance matrix will provide more accurate con- 
straints on parameter estimates. A full statistical analysis 
with marginalisation over other parameters, will be pre- 
sented elsewhere. 

An Hl-galaxy survey with a telescope such as the SKA 
provides an excellent data set for ISW and power spectrum 
measurements, both in well-defined redshifts bins. This com- 
bination may provide a unique probe of the sound speed of 
dark energy. 



ACKNOWLEDGEMENTS 

The authors would like to thank Kavilan Moodley, Bruce 
Basset and Robert Lindebaum for helpful discussions and 
comments throughout the project. ATR acknowledges the 
government of South Africa through the National Research 
Foundation (NRF) for their financial support during his 
PhD. studies and the SKA project office in South Africa 
for their support during the course of this project. 



— , 2002, PhRvD, 65, 023003 
Hu W., Jain B., 2004, PhRvD, 70, 043009 
Hu W., Scranton R., 2004, PhRvD, 70, 123002 
Lewis A., Challinor A., Lasenby A., 2000, ApJ, 538, 473 
Ma C.-P., Bertschinger E., 1995, ApJ, 455, 7 
Ma C.-P., CaldweU R. R., Bode P., Wang L., 1999, ApJL, 
521, LI 

Meyer M. J., Zwaan M. A., Webster R. L., Brown M. J. I., 

Staveley-Smith L., 2007, ApJ, 654, 702 
Nolta M. R. et al., 2004, ApJ, 608, 10 

Peebles P. J. E., 1980, The Large-Scale Structure of the 

Universe. Princeton Univ. Press, Princeton, NJ 
Percival W. J. et al., 2001, MNRAS, 327, 1297 
Perlmutter S. et al., 1999, ApJ, 517, 565 
Peroux C, McMahon R. G., Storrie-Lombardi L. J., Irwin 

M. J., 2003, MNRAS, 346, 1103 
Pogosian L., Corasaniti P.-S., Stephan-Otto C, Crittenden 

R., Nichol R., 2005, PhRvD, 72, 103519 
Riess A. G. et al., 1998, AJ, 116, 1009 
Rocha G., Trotta R., Martins C. J. A. P., Melchiorri A., 

Avelino P. P., Bean R., Viana P. T. P., 2004, MNRAS, 

352, 20 

Spergel D. N. et al., 2006, preprint ( astro-ph/0603449[ | 
Tegmark M. et al., 2004, ApJ, 606, 702 
Weller J., Lewis A. M., 2003, MNRAS, 346, 987 
Zaldarriaga M., Seljak U., 1997, PhRvD, 55, 1830 
Zlatev I., Wang L., Steinhardt P. J., 1999, PhRvL, 82, 896 
Zwaan M. A. et al., 2003, AJ, 125, 2842 
Zwaan M. A., van Dokkum P. G., Verheijen M. A. W., 
2001, Science, 293, 1800 



REFERENCES 

Abdalla F. B., Rawlings S., 2005, MNRAS, 360, 27 
Afshordi N., 2004, PhRvD, 70, 083536 
Afshordi N., Loh Y.-S., Strauss M. A., 2004, PhRvD, 69, 
083524 

Amendola L., 2000, PhRvD, 62, 043511 
Armendariz-Picon C, Mukhanov V., Steinhardt P. J., 

2000, PhRvL, 85, 4438 
Basilakos S., Kovac K., Aragon M., van de Weygaert R., 

van der Hulst J. M., 2006, preprint ( astro- ph70602073vl) 
Bean R., Dore O., 2004, PhRvD, 69, 083503 
Boughn S., Crittenden R., 2004, Nat, 427, 45 
Cabre A., Gaztanaga E., Manera M., Fosalba P., Castander 

F., 2006, MNRAS, 372, L23 
Corasaniti P.-S., Bassett B. A., Ungarelli C, Copeland 

E. J., 2003, PhRvL, 90, 091303 
Corasaniti P.-S., Giannantonio T., Melchiorri A., 2005, 

PhRvD, 71, 123521 
De La Macorra A., Melchiorri A., Serra P., Bean R., 2006, 

preprint |astro- ph/0608351 ) 
Fosalba P., Gaztanaga E., 2004, MNRAS, 350, L37 
Hannestad S., 2005, PhRvD, 71, 103519 
Hu W., 1998, ApJ, 506, 485 



